Linking confinement to spectral properties of the Dirac operator 
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We represent Polyakov loops and their correlators as spectral sums of eigenvalues and eigenmodes 
of the lattice Dirac operator. The deconfinement transition of pure gauge theory is characterized 
as a change in the response of moments of eigenvalues to varying the boundary conditions of the 
Dirac operator. We argue that the potential between static quarks is linked to spatial correlations 
of Dirac eigenvectors. 
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Introductory remarks 

Confinement and chiral symmetry breaking are two of the 
central features of QCD. At the QCD finite temperature 
transition chiral symmetry is restored and the theory de- 
confines. Numerical simulations in lattice QCD indicate 
that the critical temperature Tc is the same for both tran- 
sitions. Thus it is widely believed that there must be a 
mechanism linking the two phenomena. However, so far 
there is no generally accepted picture for such a link. 

For chiral symmetry breaking an important connection 
between the order parameter, the chiral condensate, and 
spectral properties of the Dirac operator is known. The 
Banks-Casher formula Q links the chiral condensate to 
the density of Dirac eigenvalues at the origin. Building 
on this connection it has been suggested |3] that a liquid 
of objects with topological charge, e.g., instantons, could 
lead to a non-vanishing density of eigenvalues and thus 
such topological objects are candidates for vacuum exci- 
tations responsible for chiral symmetry breaking. Above 
Tc the topological objects are expected to rearrange and 
a spectral gap opens up. The chiral condensate vanishes 
and chiral symmetry is restored. 

Concerning confinement so far no signature in spectral 
properties of the Dirac operator is known. On the other 
hand it is obvious that such signatures must exist. The 
inverse Dirac operator, i.e., the quark propagator, clearly 
knows about confinement properties. In this letter we 
present an attempt to identify spectral signatures of the 
Dirac operator which are related to confinement. 

Our starting point are Polyakov loops on a euclidean 
lattice. A Polyakov loop is defined as the ordered product 
of temporal link variables at a fixed spatial position x, 



N 



L{x) = trc J]^ Ua{x, s) , 



(1) 



i5 = l 



where the N denotes the number of lattice points in time 
direction and trc is the trace over color indices. For 
pure gauge theory the action is invariant under center 
rotations, while the Polyakov loop is not. The QCD 
phase transition can be viewed as spontaneous break- 
ing of the center symmetry and the Polyakov loop is the 
corresponding order parameter |^ . In the confined phase 
it vanishes, but above Tc, where the theory deconfines, 
the Polyakov loop acquires a non-vanishing expectation 



value. When fermions are coupled, the center symme- 
try is broken explicitly by the Dirac operator and the 
Polyakov loop cannot be used as order parameter then. 
In the dynamical case correlators of Polyakov loops are 
related to the potential between static quark sources and 
thus linked to confinement. 

Working in the lattice regularization, we express the 
Polyakov loop and its correlators as a spectral sum of 
eigenvalues and eigenvectors of the Dirac operator with 
different boundary conditions. If the Polyakov loop is av- 
eraged over space, our formula contains only moments of 
eigenvalues. The deconfining transition of QCD can then 
be seen as a change in the response of these moments to 
changing boundary conditions of the Dirac operator. We 
furthermore show, that the static potential is connected 
to spatial correlations of the eigenvectors. 

Dirac operator and Polyakov loops 

To be specific, we work on a lattice of size x N and 
for the gauge field (SU(n), n arbitrary) use boundary 
conditions which are periodic in all four directions. We 
base our discussion on the Wilson Dirac operator 



D{x\y) = {4+m)S, 
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[iTlf.] Uf,{x)6.^+ii^y , (2) 



where we use C/_^(a;) = U^{x — fiy . At the moment 
we consider boundary conditions for the Dirac operator 
which are periodic for all four directions. We remark 
that our construction goes through unchanged when a 
chemical potential is coupled in the standard way 0| . 

The hopping terms of the Wilson Dirac operator (0) 
connect nearest neighbors. When powers of D are con- 
sidered, these terms combine to chains of hops on the 
lattice. Along these chains products of the link variables 
Uf^{x) and of the matrices [1=F 7^1/2 are collected. Taking 
the m-th power will give rise to chains with a maximal 
length of m steps (in chains shorter than m the on-site 
term appears as a factor). We now consider the A'^-th 
power of D, where N is the temporal extent of our lat- 
tice. Thus we will encounter chains with a maximum 
length of N. Furthermore we set the two space-time ar- 
guments of D to the same value, y = x, such that we pick 
up only closed chains, i.e., loops starting and ending at 



2 



X. Among these are the loops where only hops in time di- 
rection occur such that they close around compact time. 
We obtain (the * denotes complex conjugation) 



trdtrc[D^(f,t|f,i)] = tr^ 



1 - 74 



N N 



trri 



1 + 74 



n N N-1 



ivc JJ^ Ui(x^N — s)'^ + other loops 



s=0 



2L{x) + 2L{x)* + other loops . 



(3) 



Here tr^ and tr^ denote the traces over Dirac and color in- 
dices. Evaluating the Dirac trace is straightforward, due 
to the projector property ([1 =F 74]/2)^ = [1 T 74]/2- We 
stress that the forward and backward running Polyakov 
loops are the only ones that wind non-trivially around 
the compact time direction. 

We now explore the fact that the Polyakov loops re- 
spond differently to a change of the boundary conditions 
compared to other, non-winding loops. We can change 
the temporal boundary condition of the Dirac operator 
by multiplying all temporal link variables at t = iV with 
some phase factor z G <C,\z\ = 1, 



Ui{x,N) -4 zUi{x,N) Vf. 



(4) 



We evaluate the left-hand side of © a second time, now 
using the Dirac operator in the field transformed accord- 
ing to which we denote as D^- We obtain: 

trd c[£'f (x,t|f,t)] = z2L(f)-)-.z*2L(x)* -I- other loops. 

. . 

Only the two Polyakov loops which wind non-trivially 
are altered when changing the boundary condition. All 
other loops cross the last temporal link, where we put 
our boundary condition, equally often in both directions 
(or not at all) . Since forward and backward oriented link 
variables are conjugate to each other the phase factors 
cancel. This fact is used to get rid of the trivial loops: 



trd.,^[D^ {x,t\x,t)\ - trrf,4L'f (f,t|f,i)] 



1 - z 



1 - z* 
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(6) 



Note that the right-hand side is independent of t and 
thus the left-hand side can be evaluated at arbitrary t. 
We use this freedom and average over t, a step which 
will be convenient later. It is easy to see that combining 
periodic boundary conditions with z = ± i boundary 
conditions (subscripts ±) gives the Polyakov loop: 

^(^) = ^ E ( 2 tr<i,c [D'^ix, t\x, t)] - (7) 
t=i 

{l+i)trd4D^{x,t\x,t)] - {l-i)trd,c[D^{x,t\x,t)]y 

We remark that the choice z = ± i is not particularly 
special and also combinations of other phases can project 



out the Polyakov loop. If one is only interested in the real 
part of the Polyakov loop, already the combination of 
periodic and anti-periodic (z — —1) temporal boundary 
conditions is sufficient as can be seen from Eq. ®. 

Spectral representation of Polyakov loops 

We now express the Dirac operator using the spectral 
theorem. Since the Wilson Dirac operator is not a normal 
operator it cannot be unitarily diagonalized and we have 
to use left- and right-eigenvectors. We denote them by vl 
and ur, with vl being a row vector, while vr is a column 
vector. They obey vl,\D — Xvl,\ and Dvr.x — Aur,a- 
The spectrum is the same for both eigenvalue problems, 
however, not the eigenvectors. The eigenvectors form a 
bi-orthonormal set, VL,i • — 5ij, and the spectral 
representation of the Dirac operator and its iV-th power 
respectively, are given by 

D = Y,XvR^x<»VL,x , = Y^X^ VR.x'»VL,x . (8) 

A A 

We remark that for a Dirac operator which is 75- 
hermitian (such as the Wilson operator we use here) the 
situation is simplified further: A 75-hermitian operator 
obeys 75-D — D^'f^, implying that the left-eigenvector 
vl.x with eigenvalue A is related to the right-eigenvector 
with eigenvalue A* (for a 75-hermitian Dirac operator 
the eigenvalues come in complex conjugate pairs) via 
vl.x — vj-^, 75, and, after suitable normalization, one 
can work with the right-eigenvectors alone. 

Inserting the spectral sum (jSJ into ^ we end up with 
our final expression for the local Polyakov loop: 

A 

il + i)Y,X^p+{x),^ - (l-i)^A^p_(f),_ 

A+ A_ 

The first sum is over the eigenvalues A obtained with pe- 
riodic boundary conditions, while the second and third 
sums are for the cases with phase factors ± i. We have 
defined the densities p and p± for these boundary condi- 
tions (color and Dirac indices are summed) as: 
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YvL,x{x, t) 
1 



VR,x{x,t) 



P±{x)x± 



t 

N 

YvL.x±(x,t) ■ VR^x±{x,t) . (10) 
t=i 

We remark that both, the eigenvalues as well as the den- 
sities p, pz are gauge invariant. Thus gauge invariance is 
explicitly manifest in 0, (|10|) as it should be, since the 
formula contains no approximations whatsoever. 

When using the Polyakov loop as an order parameter 
for confinement (pure gauge theory), it is usually aver- 
aged over the spatial sites x to improve statistics. 



P 



V3 



(11) 



3 



where V3 denotes the spatial volume. When summing 
over spatial indices as well, the densities H1(J|) turn into 
the matrix elements of left- and right-eigenvectors with 
i = j and thus the summed densities equal 1. We obtain 
for the averaged Polyakov loop, 

^ - ^ f 2 E - (1 +i) E - (1 -i) E V 

\ A A+ ^- / 

(12) 

where = L'^ is the total number of lattice points. 
The equation expresses the Polyakov loop through mo- 
ments of the Dirac eigenvalues at different boundary con- 
ditions. We remark that the individual sums are real 
since the eigenvalues come in complex conjugate pairs. 

Let us add a few remarks on the case when an exactly 
chiral lattice Dirac operator is used. A chiral operator 
obeys the Ginsparg Wilson Isl equation, 75Z) -I- D75 = 
D"f^D, which governs chiral symmetry on the lattice. All 
the properties of the Dirac operator we have used so far 
still hold. There is, however, an important qualitative 
difference when using a Ginsparg Wilson Dirac opera- 
tor. A solution of the Ginsparg Wilson equation is not 
ultra-local, i.e., D{x\y) contains not only nearest neigh- 
bor hops, but paths of all lengths. The coefficients for the 
paths decrease exponentially with their length. The exis- 
tence of paths of all lengths implies that will not only 
give rise to non-trivially winding loops that are straight 
lines in time direction, but also "dressed" Polyakov loops 
taking some detour in spatial directions. These dressed 
loops are exponentially suppressed with length (as long as 
one keeps away from the Aoki phase - see [g)- However, 
also the dressed Polyakov loops are not invariant under 
Zjv transformations and thus are equally well suited or- 
der parameters as are the straight Polyakov loops. 

Discussion of the spectral sums 

We begin the interpretation of our formulae with the 
spectral sum (|12f) for the averaged Polyakov loop P. 
The right-hand side expresses the Polyakov loop as lin- 
ear combination of the A'^-th moments of the eigenvalues 
of the Dirac operator at different boundary conditions. 
When crossing the QCD phase transition the Polyakov 
loop acquires a non-vanishing expectation value (pure 
gauge theory), and the interplay of the eigenvalues at 
different boundary conditions must change at Tc- 

An interesting question is which part of the Dirac spec- 
trum carries most of the signal for the Polyakov loop. 
The infrared part of the spectrum (small eigenvalues) is 
known to undergo a pronounced change as one crosses 
from the confining to the deconfined phase: A gap opens 
up in the spectrum, the density of eigenvalues near the 
origin vanishes and, according to the Banks-Casher for- 
mula, chiral symmetry is restored. The size of the gap 
depends on the relative phase of the Polyakov loop and 
the Dirac boundary condition . 

Although the low-lying eigenvalues undergo a dramatic 
change, it is not clear whether they give a sizable contri- 



bution to the Polyakov loop: Due to the large power 
the small eigenvalues are strongly suppressed relative to 
the bulk of the spectrum where eigenvalues of 0(1) occur 
(in lattice units). On the other hand, since the hopping 
term changes sign under a staggered sign transformation, 
for each eigenvalue A also 8 — A is an eigenvalue. These 
mirror images undergo the same drastic change at the 
phase transition, but are not suppressed. Thus small 
eigenvalues could contribute indirectly through their mir- 
ror images at the cutoff. A second important aspect is 
that the density of eigenvalues increases roughly as |Ap. 
Thus the infrared part of the spectrum and its mirror 
image at 8 — A have only a small weight. How much each 
part of the spectrum contributes to the spectral sum ifT^ 
can probably only be decided numerically. 

In a first test we summed up the lowest 100 eigenval- 
ues of the chirally improved lattice Dirac operator for 
quenched SU(3) configurations on 20^ x 6 lattices. This 
was done for temperatures below and above Tc and in 
both cases summing up the low-lying modes gave only 
a small fraction of the true Polyakov loop P as deter- 
mined directly from the link variables. It is remarkable, 
however, that in all cases we checked for T > the 
spectral sum (|12|l evaluated with even only the 50 lowest 
eigenvalues already gets the phase of the Polyakov loop 
right (the Z3 symmetry gives rise to 3 preferred phases 
of the Polyakov loop expectation value) . Unfortunately 
very little is known about the ultraviolet part of the spec- 
trum. Numerical evaluation of the eigenvalue spectrum 
is usually restricted to the infrared part or to very small 
lattices due to the immense computational cost of a com- 
plete diagonalization. Random matrix theory, which is 
very successful in describing the infrared part of the spec- 
trum, is valid only up to the Thouless energy and pro- 
vides no insight to the ultraviolet part of the spectrum. 

We now turn to interpreting Eq. ^ which is the spec- 
tral sum for the local Polyakov loop at x. In its local 
form the Polyakov loop L{x) can be used for computing 
the potential V{r) between static sources (r = \x — y\), 

{L{x) L{yY) oc exp ( - iV V{r)) - exp ( - aNr) , (13) 

where the last expression is the behavior for large r, with 
a denoting the string tension. When inserting the spec- 
tral sum lO in the correlator l|13|l . one sees that the 
static potential is related to correlations of the densi- 
ties p and p±. In particular correlations of the form, 
{p{x)\p{v)y). {p{x)xp±iy)yj, {p±{x)x^p±{y)*yj, build 

up exTp{— N V {\x — y\)) . In the confining phase the string 
tension a is non- vanishing, such that the density corre- 
lations are expected to decay exponentially with |a; — y| 
(until string breaking sets in). In the deconfined phase 
one has a = and the decay of the eigenvector density 
correlators should show a power law behavior. 

Finally, let us address another possible application of 
our result: The spectral representation (O of L{x) might 
be an interesting observable also for individual gauge con- 
figurations. It provides a filter for analyzing topological 
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infrared structures of the gauge field. Expressing glu- 
onic observables in terms of spectral sums and truncat- 
ing these sums has been applied before (for examples see, 
e.g., The technique is sometimes referred to as "low 
eigenmode filtering". It makes use of the fact that the 
low-lying modes of the Dirac operator are an efficient 
filter for infrared properties of the gauge field. In a ther- 
malized configuration, i.e., one not treated with cooling 
or smearing, these are hidden under UV fluctuations and 
a filter is needed to observe them. Using the low-lying 
Dirac modes has the advantage over smearing or cooling 
techniques, that the gauge field is not altered. Also low- 
lying modes of the covariant Laplace operator have been 
implemented as a low pass filter 10]. 

The newly developed formula for the Polyakov loop 
might be particularly useful for analyzing properties of 
so-called Kraan-van Baal solutions [ll|. There an object 
of topological charge 1 is made from n constituents (for 
gauge group SU(n)) and the local Polyakov loop is ex- 
pected to show a specific pattern at the positions of the 
constituents. For cooled configurations this pattern has 
been confirmed j^. The expression ^ can be used to 
study the local behavior of the Polyakov loop directly for 
thermalized configurations using low eigenmode filtering 
techniques. Concerning the dominant contributions, for 
the local version L{x) the situation is different from the 
summed Polyakov loop P. In the spectral sum ^ the 
eigenvalues are multiplied with the densities p{x)x and 
p{^)\± ■ For the low-lying eigenvectors these densities are 
known to be large for positions where topological objects 
are located, while for larger values of A, \± the densities 
are dominated by fluctuations. Numerical tests of using 
@ as an infrared fllter are in preparation. 

Concluding remarks 

How chiral symmetry breaking reflects itself in spectral 
properties of the Dirac operator is well understood. The 



chiral condensate is related to the density of eigenvalues 
near the origin and random matrix theory can be used 
to describe the behavior of the low- lying eigenvalues. 

Also confinement should have a signature in spectral 
properties of the Dirac operator, which, however, is only 
vaguely understood so far. In this letter we have shown 
that the Polyakov loop can be expressed through mo- 
ments of the Dirac eigenvalues computed with different 
boundary conditions. At the phase transition of pure 
gauge theory the Polyakov loop acquires a non- vanishing 
expectation value and the dependence of these moments 
on the boundary condition changes. We also discuss 
correlators of local Polyakov loops which are related to 
the potential between static quarks. Based on this re- 
lation we suggest that in the confining phase correlators 
of Dirac eigenvector densities decay exponentially with 
spatial separation, while a power law should be seen in 
the deconfined phase. 

Some of our spectral sums might get sizable contri- 
butions from large eigenvalues. Analyzing the ultraviolet 
part of the Dirac spectrum is a challenge both for numeri- 
cal and analytical approaches. The results presented here 
are certainly only an initial step indicating which spec- 
tral quantities can be studied. The ambitious goal of such 
a study could be a link between the spectral properties 
in the infrared, which are connected to chiral symmetry 
breaking, and the ultraviolet part probably more relevant 
for confinement. 
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